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Abstract 

We determine non-perturbatively the parameters of the lattice HQET Lagrangian and 
those of the time component of the heavy-light axial-vector current in the quenched 
approximation. The HQET expansion includes terms of order 1/mt,. Our results allow 
to compute, for example, the heavy-light spectrum and B-meson decay constants in the 
static approximation and to order in HQET. The determination of the parameters 

is separated into universal and non- universal parts. The universal results can be used 
to determine the parameters for various discretizations. The computation reported in 
this paper uses the plaquette gauge action and the "HYP 1/2" action for the b-quark 
described by HQET. The parameters of the current also depend on the light-quark 
action, for which we choose non-perturbatively 0(o)-improved Wilson fermions. 
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1 Introduction 



Heavy quark effective theory (HQET) was developed already quite a while ago [T]-[5] . 
Still it is of considerable interest today, primarily for two reasons. First it describes 
the asymptotic expansion of QCD observables in the limit of a large quark mass, in 
particular the mass of the b-quark, m^. For this to be true, the observables have to be in 
the proper kinematical region. But in such a region, we expect the expansion in 1 jm\, to 
be valid also non-perturbatively if the parameters in the effective theory are determined 
non-perturbatively. Understanding QCD then includes understanding the HQET limit. 
Steps towards establishing the expected equivalence on the non-perturbative levej]]have 
already been carried out |6,7|, but we intend to go much further, in particular through 
a complete treatment of 1/rrib corrections. 

Second, flavour physics has become a precision field. No flavour physics observ- 
able has shown evidence for physics beyond the Standard Model with the presently 
available precision. In many cases the limiting part is the theoretical uncertainty, not 
the experimental one [:8|. Reliable lattice computations are needed to make progress, 
but heavy quarks are difficult due to 0((amb)") discretization errors, where a is the 
lattice spacing. In particular, if am^ is too large the expansion in a will break down 
altogether [9,10 Instead, in HQET the cutoff effects are 0((aAQCD) n ); they are thus 



much more manageable. 



matching SSF 
QCD(Li) " HQET(Li) HQET(L 2 ) 



Figure 1: The strategy. 

Here we present the determination of the parameters of the effective theory in 
a formulation where all power divergences are subtracted non-perturbatively through 



a matching of QCD and HQET in a small volume 16 with Schrodinger functional 
boundary conditions. As illustrated in Fig. [TJ we then continue to larger lattice spa- 
cings through a step scaling method. Section [3] explains the details, in particular how 
lattice spacing errors are removed in each step. The work of JTt] is here extended by a 
determination of all parameters in the action as well as those for the weak current Aq 
including the terms of order l/m^. 



1 Note that to our knowledge a proof of the renormalizability of HQET to all orders of perturbation 
theory has not yet been given. 

2 The second one of these references shows the difficulty even for a charm quark. For more general 



discussions and more references we refer the reader to reviews 11 15 
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Our numerical implementation is done in the quenched approximation. It repre- 
sents a very useful test-laboratory, in particular for the first motivation presented above: 
the qualitative features of the 1 /m^ expansion will not depend on the fact that the light 
quarks are quenched. Applications of the parameters computed in this paper require 
the determination of energies and matrix elements in a large volume and will appear 
in separate papers, but already here we learn interesting lessons about the asymptotic 
convergence of the expansion. In the following section we recall the basic formulation 
of HQET, with a primary focus on defining the parameters of the theory and how they 
enter the computation of the spectrum and matrix elements. Section 3 gives a short but 
precise account of the strategy for the computation of the parameters. The technical 
Sections 4.1-4.2 discuss the numerical details of the different intermediate steps while 
Section 4.3 presents our main results for the HQET parameters. 



2 HQET at order l/m b 

In this section we define HQET including terms of order l/mb, m particular the pa- 
rameters of the theory. We then show the expansion of a few observables as examples 
how the parameters can be used. We choose to regularise the theory on the lattice 
although our approach is in principle independent of the specific regularisation. Almost 
all formulae are established in [17] . They are repeated here for the benefit of the reader 
to keep the paper self-contained, but for details, such as the exact consequence of spin 



symmetry, the reader is referred to 17 . Terms of order 1/m? are dropped without 
notice. 



2.1 Lagrangian 

The HQET Lagrangian, 

•^HQETOr) = £?sta,t(x) - UJ khl O kin (x) - LJ spin O spin (x) , (2.1) 

consists of the lowest order (static) term, 

^stat(x) = i> h {x) DoMx) , (2.2) 
and the first order corrections 

O ki n(x) = ^ h (x)D 2 V h (x) , O spin {x) = ^ h (x)a-B^ h (x). (2.3) 



We use the backward covariant derivative Dq as in [18] and the 4-component heavy 
quark field subject to the constraints P+iph = iph , ip^P+ = V'h with P + = (1 + 7o)/2 
and the discretized version ctB = ^ k ■ akjF^j / (2i) , where Okj and the lattice (clover) 
field tensor F are defined in 19 . The kinetic term D 2 is represented by the nearest 



neighbour covariant 3-d Laplacian. The normalization is such that the classical values 
of the coefficients are u: k - m = uj spm = l/(2mt>). A bare mass m-bare has to be added to 
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the energy levels (e.g. the B-meson mass) computed with this Lagrangian to obtain the 
QCD ones (up to 1/m^). At the classical level it is m^, but in the quantized theory, it 
has to further compensate a power divergence. 

2.2 Weak axial current 

For many applications, the weak, left handed, heavy-light current is needed. We here 
consider just the time component of the axial-vector part, Aq. The other components 
can be treated analogously, but so far we restricted ourselves to Aq, which is sufficient 
for a computation of pseudoscalar decay constants. 

At the lowest order the current is form-identical to the relativistic one. At first 
order it is corrected by two composite fields of dimension four. The explicit form is 

A^(x) = zf QET [Af*\x) + J24 ) A$\x)] , (2.4) 

i=l 

A$\x) = ^(x)^757i(Vf-^^ h (x), (2.5) 

4 2) (x) = -diAf^ix), ^f at (x)=^ 1 (x) 7i7 5V'h(x), (2.6) 

where all derivatives are taken to be the symmetric nearest neighbor ones, 

d i = \(d i + &t), 5f= + Vf=|(Vi + Vt). (2.7) 

By considering Z^ ET to be a function of m\/m\, with the light quark masses mi, 
the above set of operators is complete after using the symmetries and the equations of 
motion. Since m\ and such effects are further reduced by a factor of the coupling 

constant a{m^) we ignore this dependence and determine Z^ ET with the light quark 

(2) 

mass set to zero. Note that Aq (x) does not contribute to correlation functions and 
matrix elements at vanishing space-momentum. It is not needed for the computation 



of decay constants and has not been written down in 17 . We will also not determine 



its coefficient here. 

A short note on discretization errors is useful at this point. Remaining in the static 
approximation (wkin = ^spm = 0), the Lagrangian has been shown to be automatically 
0(a) improved by studying its Symanzik expansion. Therefore, all energy levels are 
as well Il9|. The zero space-momentum matrix elements of the current A^ ET are 
0(a) improved in static approximation if one sets = ac^ at 18 19j. |^] Including the 
1/mb terms, linear terms in a remain absent, except for those accompanied by a factor 
1/mb 16 . So leading discretization errors are 0(a/mb) and 0(a 2 ). 



3 At zero space momentum, partial summation can be used to bring Aq 2 ' into the form used in 19 
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2.3 Observables 

2.3.1 Correlation functions 

Observables of interest are obtained from Euclidean correlation functions in large vo- 
lume. As an illustration we consider the QCD correlator 



Caa(xo) = Z 2 A a 3 ]T (a (X) (A (0)) 



of the heavy-light axial current in QCD, = "017^75^ (and = — ^7^75^1)- ^ * s 
normalized by Z\ to satisfy the chiral Ward identities 20 2l]. Ignoring renormalization 
for a moment, an expansion in reads 

Caa(x ) = Ctl\x ) + c£l m \x ) = Ctl\xo)[l + R i li m \xo)} (2.9) 

where the notation indicates Caa = 0(V m b)- I n HQET the path integral weight is 
expanded in l/m^ and one obtains the fully renormalized expansion 

Caa(x ) = e - m ^o (z HQET )2 [ c stat (xo) + UkinC fl( Xo ) + a; spin C A p A n (x ) 

+4 1) [Qaa(^o) + Cff A (x )]] (2.10) 

1 + a; kin i?^(xo) + io spin R s l A n {x ) + c A X) ^(xq) 
in terms of bare static expectation values, for example 

C^(xo) = a 7 ^(< at (x)« at (0)) t Okin(^)) • (2.12) 

' \ / stat 

X, z 

2.3.2 Spectrum and matrix elements 

The spectral representation yields the large time behaviour in QCD 

Caa(xo) = A 2 e~™ [ 1 + 0(e" A * ) ] , (2.13) 

where A is a gap in this channel ( ~ 2m n in large volume) and A = (B(p = 0)|^4o(0)|0) 
in the non-relativistic normalization of states {B(p)\B(p')) = 2S(p — p') . Of course 
A = /b y /?ti b/2 is a phenomenologically interesting quantity. Again a naive expansion 
reads 

CaaOco) = (A stat ) 2 e-< atxo [l + 2A (1/m) -m£ /m) x ]. (2.14) 

Using the transfer matrix of the static theory, the HQET correlators (static and beyond) 
are easily analysed. The result is 

C$£(x ) = (A stat ) 2 e - EBt * txo [ 1 + 0(e" Astata;o ) ] (2.15) 
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for the static term and for example 

Rfl(x ) = 2A kin - x E kiQ + O(x e" Astat x °) (2.16) 
for the kinetic correction^ 



A comparison to eq. ( |2.13 | yields 



m B = m hare + E st ^ + cu kin E kin + u spin E spin , (2.18) 

ln(A rl /2 ) = ln(2 A IQET ) + ln(„4 stat r 3 Q /2 ) + A 5A + w kin ^ kin + a; spin ^ spin , 

(2.19) 

where ro is an arbitrary length-scale of the theory. Completely analogous formulae hold 
for the excited state masses and matrix elements. The necessary bare quantities such 



as E km can be efficiently computed with the generalized eigenvalue method 22 . We 



o in 

have here chosen to write the expansion of In (-4 r ) instead of A itself, since in this 



way one explicitly avoids terms quadratic in 1/mb, while e.g. in eq. (2.10) terms such 
as Z^ mh '' 'wklnCAA are understood to be dropped to remain consistently in order 1/mb- 
This rule is necessary for a correct renormalization of the theory. 



2.4 Properties of non-perturbative HQET parameters 

Before entering the presentation of a computation of the HQET parameters, we here 
mention some of their properties. 

2-4-1 Scheme dependence 

Since we are working in the framework of an effective field theory, which has non-trivial 
renormalization, the parameters have to be determined by a matching performed at 
finite value of the expansion parameter 1/mb- The truncation of the effective theory 
then introduces an ambiguity of order of the left out terms. 

As a result the parameters such as u>km an d Z A ^ ET have a dependence on the 
choice of the matching condition. This is a scheme dependence, analogous to the one 
of the usual perturbative expansion in the renormalized QCD coupling. For an explicit 
example consider ln(Z A ). Evaluated in the static approximation we will denote it 
by ln(Z A tat ) and including all l/m b terms it is denoted as ln(Z A QET ) = ln(Z A tat ) + 
ln(Z A 1//mb ^). The scheme dependence of ln(Z A ^ ET ) is then of order l/?n?, while the 
one of ln(Z A at ) and ln(Z A ^ mb ' ) ) individually is of order 1/mb- Note that in all this 
discussion we are working non-perturbatively in the QCD coupling but order by order 
in 1/mb- 



4 An explicit expression for E kln is 

E kin = -^(B|a 3 ^O kin (0, Z )|B) stat =-i<B|Oida(0)|B) rtrtl (2.17) 

z 

where the state \B) is an eigenstate of the static transfer matrix. 
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2-4-2 Difference to previous renormalizations in the static approximation 



In order to avoid a misunderstanding we also point out the difference to the renormal- 



ization carried out for example in 23,24 . As long as one is working just in the static 
approximation, the weak currents do not mix with operators of different dimensions. 
It is then also consistent to renormalize them perturbatively. Both to avoid ambigui- 
ties in renormalization schemes and to be able to profit from the high order continuum 
perturbative results (25l[26) , it is advantageous to first introduce the RGI current, 



(A* Gl ) (x) = Zt% G1 Ar t (x). (2.20) 
The renormalization constant has been determined non-perturbatively in 



23 



for the quenched case, and in 27 for two flavours of dynamical quarks. 

However, here we determine Z^ at by a non-perturbative matching to QCD, and it 
is different. In fact the correspondence is 

Z stat = z$$ G1 C PS (M h /A) , (2.21) 
in terms of the matching function Cps introduced in [24] . In the "old" strategy of for 



example 24, 28 , functions such as Cps are determined from high-order perturbation 
theory, while here we evaluate the full factor Z^^ non-perturbatively. 

3 Computation of HQET parameters 





definition 


classical 
value 


static 
value 


U\ 


^bare 


m b 


m bare 


Ldo 


ln(zl QET ) 







L0 3 


c (1) 


-l/(2m b ) 




LOl 


<^kin 


l/(2m b ) 





tor, 


^spin 


l/(2m b ) 






Table 1: Notation for HQET parameters, their values in classical and static approximation. 

We collect all HQET parameters into one vector u with components Ui , % = 1, . . . , 5 
listed explicitly in Table [TJ In static approximation the parameters are u>i, u)%- When 
l/m b corrections are included, the additional parameters in the action are W4, W5; more- 
over the previous parameters change by (partially power divergent) terms of order l/m b . 
The situation with W3 is more intricate. It is needed for 0(a)-improvement of the static 
approximation and for genuine l/m b -terms at order l/m b . We now turn to an expla- 
nation of the various steps involved in the determination of the Wj. Our strategy is 
illustrated in Fig. [2} 
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HQET 




11 



35 



5' 5 



Figure 2: Illustration of the strategy. Our numerical application uses much finer resolutions a/L than 
those shown here. For each step i a series of simulations Si is necessary, they are described in the text 
and in Table 1] 



3.1 Finite volume observables 

Five observables $j are required to determine the HQET parameters aj{. We choose 
&i(L, M, a) defined from Schrodinger functional correlation functions by forming suit- 
able renormalized combinations. They are universal, which means in particular that 
their continuum limit &i(L, M, 0) exists. 

As variables we have chosen the box size L (which plays the role of a kinematical 
variable), the RGI mass of the heavy quark, M, and the lattice spacing a. Equivalent 
but dimensionless variables are a non-perturbative running coupling g(L) (e.g. in the 
Schrodinger functional scheme), the combination z = ML and the resolution a/L. The 
light quark masses are assumed fixed (in the numerics they are set to zero). 

Most of the details of our choice follows closely |17| ; for others we refer the reader 
to App. [A} Here we note only a few properties. As L becomes large, , i = 1, 2 tend to 
the B-meson mass and the logarithm of the B-meson decay constant, up to kinematical 
constants. They are thus mainly determining Wj<2- Correspondingly $4 and $5 allow 
easy access to the 1 / mb-parameters in the Lagrangian and <3?3 to the correction term of 
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the current with coefficient cjj . 

By a choice of notation (e.g. considering ln(„4.) instead of A), we arranged our finite 
volume observables to be linear in cj, 

M, a) = 7](L, a) + <p(L, a) u(M, a) . (3.1) 

By construction <f> is a block diagonal matrix, which is explicitly given in terms 
of (bare) correlation functions in HQET. The inhomogeneous pieces iji involve just 
correlation functions in the static approximation. Their continuum limits rji(L,0) exist 
for % > 2, while for i = 1,2 additive renormalizations are necessary. They are contained 



in the second term of eq. (3.1). In the case i > 2, the two terms of the right hand side 



of eq. (3.1) are then computed separately, see section 3.5 

The HQET parameters oj are defined by matching at a certain value of L = L\ the 
HQET observables <3?j to the corresponding values $9 CD computed in QCD. These QCD 
observables are first extrapolated to the continuum limit (indicated by S\ in Fig. [2J 

$, QCD (Li, M, 0) = lim $? CD (Li, M, a) . (3.2) 

i a— >0 

3.2 Matching 

The matching scale L\ must be chosen such that 1/L\ ^ to allow for a precise 
expansion in l/mb- Moreover, we want lattice spacings of order 10 -2 fm in order to keep 



omb < 1/2 while performing the continuum extrapolation eq. (3.2). These constraints 



lead to a box size L\ ~ 0.4 fm 17,24 



The HQET parameters are then defined by imposing 

$(£i,M,a) = $ QCD (Li,Af,0) (3.3) 



for any value of the lattice spacing. By inverting eq. (3.1) one obtains 

oo(M,a) = 0- 1 (Li,a)[$ QCD (L 1 ,M,O)-r ? (L 1 ,a)]. (3.4) 
Through our choice the matrix (j) is of the form 



C B\ j _ fC- 1 -C- l BA~ l 

A ' " V A- 1 



(3.5) 



where C = diag(L, 1). The 3x3 sub-matrix A has a further upper triangular structure. 
It is written explicitly in App. |Bj together with its inverse. 

3.3 Step scaling 

For reasonable resolutions L\/a > O(10), only lattice spacings a < 0.05 fm are accessi- 
ble, while for standard large volume HQET computations we would also like larger a. 
In order to obtain u also there, we first need 3>(L, M, 0) at larger L = L2. We simply 
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use the values u(M,a) of eq. (3.4) and determine the continuum limit of the HQET 



observables at L 2 (indicated by S3 in Fig. [2]) 

$(L 2 , M, 0) = lim { V (L 2 , a) + 0(L 2 , a) w(Af, a)} . (3.6) 

a->0 

This can be done as long as the lattice spacing is common to the n 2 = L 2 /a and 
n\ = Li/a-lattices and 

s = L 2 /L\ = n 2 jn\ (3.7) 

is kept at a fixed, small, ratioj^] We choose s = 2 in the numerical application. Typical 
resolutions are L\/a = O(10), see Sect. [4j This procedure which takes us from L to 
sL with a finite scale factor s is called step scaling [ 29] . The explicit forms of the step 
scaling functions are given in App. [Bj 

3.4 HQET parameters 



The parameters Ui for use in large volume (see Sect. |2.3.2) are finally obtained from 



u(M, a) = (f>-\L 2 , a) [$(L 2 , M, 0) - rj(L 2 , a)] (3.8) 

through HQET simulations (indicated by S4 in Fig. [2]) with typical resolutions L 2 /a = 
O(10). This results in lattice spacings of a ~ 0.1 fm and smaller, as desired for the ap- 
plications. While the final results cu(M, a) do depend on all details of the regularization 
(gauge action etc.), the intermediate observables <&i(L 2 , M, 0) may be used universally 



(e.g. together with an implementation of eq. (3.8) with a different regularization). 



3.5 Splitting lowest order and first order in l/m\, 



In Sect. 2.4. 1| we discussed that the splitting of a prediction into different orders of the 



1/mb-expansion is not unique. Nevertheless, it is of interest to organize the calculation 
into a static one and the remainder. In this way one can judge on the generic size 
of 1/m-b corrections and thus get some indication of the asymptotic convergence of 
the series. A second reason is that the static theory is O(a) improved when U3 is 
included and then represents the improvement coefficient ac^ at = 0J3 |l9| . Thus the 
by far dominating part of the result can be extrapolated quadratically in a and only a 
small correction has to be extrapolated linearly to the continuum limit. Since a non- 
perturbative determination of c^ at has not been carried out, we will here use its one- loop 



perturbative value 30 . We anyway carried out quadratic extrapolations in a and then 
studied the effect of incomplete improvement. We will see in Sect. |i| that cf at is of 
little relevance. Therefore we profit from the 0(a)-improvement of the static theory, 
irrespective of the precise value of c^ at . 

5 The ratio s is required to be fixed for the cutoff effects to be a smooth function of a/ Li. 
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To summarize, the static approximation is denned by using exactly the same for- 
mulae as in HQET with l/m^ but setting 



w f at = acf at (3.9) 

^stat = w stat =Q; $|tft = Q ^ 310 ) 

All quantities are determined once in this static approximation and once with 1/m-b 
terms included and then the pure 1/mb-corrections are given by 

$(l/m) _ $ _ $ stat ? ^(l/m) _ ^ _ ^stat _ (3.H) 

For any new observables, in particular in large volume, the static results are given 



by inserting o; stat into the expressions such as eq. (2.19) and the 1/mb-correction by 
inserting u^ 1 ' m ' instead. Since, as we discussed earlier everything is to be linearized in 
the co, the "full" result up to l/?rt2 is obtained from summing static and 1/mb-correction 
or by using directly uj. We finally note that even though the "full" result contains 0(a) 
discretization errors, it appears justified to extrapolate numerical data with a leading 
correction term oc a 2 since the linear terms are suppressed by a small 1 /my, factor, which 
one would estimate e.g. as l/(mb?"o) ~ 1/10. 

4 Numerical results 

The numerical computations have been performed in the quenched approximation. For 



the QCD part we used non-perturbatively 0(a)-improved Wilson fermions [31], see 17 



for details of the Schrodinger functional implementation. Two discretizations of the 



static quarks have been considered: the so-called HYP1 and HYP2 actions [18 32 . In 
order to reduce discretization effects, we have implemented tree level improvement of 
our QCD observables and of the HQET step scaling functions (see App. [Pj ). 

Our observables depend on three periodicity angles Oq,9i,02, (see App. [A]). As 
in 17 we considered all combinations of those with 9i 6 {0, 0.5, 1.0} , 6\ < 62, i.e. nine 
different matching conditions. For our discussion of the results and the extraction of 
the final HQET parameters we chose a "standard set" 

= 0.5, 6»i = 0.5, 6» 2 = 1.0 (4.1) 

because this yields overall the smallest statistical errors in the HQET parameters. We 
will comment on the spread of results with other choices for 9i as we go along. 

We summarise the different simulations needed in our strategy in Table [2] Since a 



size L = L\ is used for the matching (eq. (3.2), eq. (3.4)) both QCD (Si) and HQET 



(S2) are simulated in that volume. The volume of space extent L2 was simulated (in 
HQET), with two different sets of lattice spacings: 



First (53) with the same set of lattice spacing as used in S2, to evaluate eq. (3.6). 
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Simulation 


L 


Theory 


L/a 


Si 


Li 


QCD 


40,32,24,20 


5*2 


Li 


HQET 


16,12,10,8,6 


s 3 


L 2 


HQET 


32,24,20,16,12 


S4 


L 2 


HQET 


16,12,8 



Table 2: Summary of the simulations used in this work. Note that L\/a — 16, L-zja = 32 are in 



addition to those of 17 



Then (S4) with lattice spacing of order 0.05 ... 0.1 fm, in order to compute u(M, a) 



given by eq. (3.8) 



The choice of the simulation parameters is described in detail in |17| ; Table 3 of that 
paper lists those for Si, Table A.l of [16] the parameters of S2 and S3 and Table 6 
of [18] the parameters of 5*4. Here we just note that L2 = 2L\ and L\ is fixed by 
the Schrodinger functional coupling g 2 (L\) = u\ = 3.48. Since this condition was 
implemented only within a certain precision, a small mismatch of g 2 (L) = u\ used in 
5*1 and g 2 (L) = ui used everywhere else has to be taken into account. This is done in 
complete analogy to appendix D of [17]. Without discussing the details we will quote 
the small corrections below. 

4.1 Observables in L = L\. 

The considered observables naturally split into three classes according to their magni- 
tude in the 1/mb expansion, or more precisely their order in z~ l = (ML) -1 . We have 
$1 = O(z), $2. ..4 = 0(1) and $5 = 0(1 /z). As an illustration of the numerical determi- 
nations, in particular the continuum limit, we discuss examples of all three cases. The 
precise definitions of $i are found in appendix |A{ 

We recall that $1 is simply a finite volume pseudo-scalar meson mass (up to 
a normalization) and determines the quark mass. In Fig. § we show its contin- 
uum extrapolation for our chosen values of the heavy quark mass, corresponding to 
z = L\M = 10.4, 12.1, 13.3. We have taken into account the errors coming from the 
relation between the bare quark mass and the RGI quark mass. In particular, a part 
of this error is common to the data at all lattice spacings. It is included after the 
continuum extrapolation; the final error bar is shown on the left side of this plot (see 
also |17|). Since O(a) improvement is implemented, this as well as all other extrapola- 



tions of eq. ( pT2| ) are done with the form ®i(Li,M,a) = $i(Li,M,0) + Cj(M) (a/Li) 2 . 

Out of the O(z ) observables, $2 determines the normalization of the axial current. 
Its continuum extrapolation, eq. (3.2), is also illustrated in Fig. [3]. 



Some of the static quantities rji (L\ , a) have a well defined continuum limit (see 



appendix [B[). In these cases, rji(Li,a) is replaced by r/j(Li,0) in eq. (3.4). We show 



examples of the determination of the continuum limit rjs(Li, 0) and 774 (L±, 0) in Fig. [4j 
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Figure 3: Continuum extrapolation of <J > ^ CD (Li, M, a) and <3?^ CD (Li, M, a), for the three different 
quark masses (from the lightest to the heaviest: upward-pointing, right-pointing and downward-pointing 
triangle). As explained in the text, the final errors for $? CD (Li,M, 0) are shown on the very left. Only 
the data in the range of the full lines are used in the continuum extrapolations. 



The graphs use the one-loop values of the improvement coefficient 30 



stat 



stat 



0.0029 gl for the HYP1 action . 
0.0518 gl for the HYP 2 action , 



(4.2) 



As a test in how far the remaining uncertainty in c^ at matters, we repeated the calcu- 
lation setting c^ at = instead. The change is invisible in the graphs and in general, 
for all the continuum extrapolations of this paper, the changes in continuum values are 



much smaller than our quoted errors. We thus consider eq. (4.2) entirely sufficient 
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Figure 4: Continuum extrapolation of 773 (Li, a) = i?A at (I/i, a), and l]i{Li, a) = Rf at (Li, a). The three 
different choices (^1,^2) = (0,0.5), (0.5,1), (0,1.0) are represented by squares, circles, and diamonds 
respectively. In each case two different static actions, HYP1 (the red points, slightly shifted to the right 
for visibility) and HYP2 (the blue points) are used. The continuum limit is obtained by a constrained 
fit rji(Li,a) — r)i(Li,0) + a,j a 2 /L\, with j = 1,2 for the two different actions. 



We turn now to $4 and $5 , needed for the determination of the 1 / my, parameters 
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Wkin and w sp i n . Their continuum extrapolations are shown in Fig. [5] Due to the exact 
spin symmetry of the static effective theory, $5 has no static contribution. Note that 
both $5 and the pure l/m^ part of $4, obtained after subtraction of 774(^1,0) (see 
Fig. [4]), are an order of magnitude smaller than <I>4. As expected $4 — 774 and <3?5 are 
decreasing functions of the quark mass. 
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Figure 5: Continuum extrapolation of the QCD observables $J CD (_Li,M, a) and $^ CD (iji, M, a) for 
the standard 9's. The three different quark masses are shown, with the same conventions as in Fig. [3] 



4.2 Observables for L = L2. 



The step to observables in the larger volume is described by eq. (3.6). It can be broken 
up into several step scaling functions defined in appendix [Bj which individually have a 
continuum limit. Discussing them one by one would be too lengthy and is also not very 
illuminating. We follow exactly Sect. [3] and insert eq. (3.4) into eq. (3.6) lattice spacing 
by lattice spacing and then extrapolate a I Li — > 0. 

Two examples are shown in Fig. [6j We observe that having the data for two 
different static actions is very useful to constrain the continuum limit, particularly so 
for the 1/mb parts. Also the resolution Lija = 32 which is in addition to those of 17 
helps a lot (the reader is invited to compare $^ 1//m ^ in Fig. [(3] to Fig. 5 of 17 ). In all 
cases the 1 /rrtb corrections are much smaller than the leading terms, suggesting a good 
asymptotic convergence of the 1/mb expansion. The precision of the 1/mb term is worse 
than the static one for $2, partly since the latter can be extrapolated quadratically in 
a to the continuum. For $1 this is different because the overall error contains a large 
piece coming from the renormalization factor determining the RGI quark mass in QCD. 

For reasons of numerical precision, $5 (and only $5) is not computed exactly as de- 
scribed in Sect.[3j Its definition (App.|A]) involves the propagation of a heavy quark over 
a distance T = L, introducing significant statistical errors for large L/a in the effective 



theory. These become unpleasantly large in eq. (3.6), more precisely in ^55(^2, a). We 
therefore replace $5(1/2) by (^(Z^) differing only by the choice T = L/2. An obvious 
question is why this is not done already for L = L±. The reason is that $5(^1) turns 



13 



18.5 

d 18 

+= 

01 1 — I 

17.5 





0.8 




0.6 




0.4 




0.2 


T— 1 




> — ' 









-0.2 




-0.4 




0.02 0.04 

a/L 2 



x 10 



4 

-3 



0.06 



0.6r 

0.58 
0.56 
& 0.54 
0.52 
0.5 



-f-j 

02 <M 





0.06 




0.04 


? 






0.02 


- eg 















-0.02 



2 



0.02 0.04 

a/L 2 



x 10 




0.06 



Figure 6: Continuum extrapolation of <J?i(Z/2) and $2^2) for the central mass z = L\M = 12.1 and 
for the standard 8's. The scale of the y-axis of static (top) and (1/m) parts is chosen to be equal. As 
explained in the text, the static parts are extrapolated quadratically in a while the (1/m) contributions 
are extrapolated linearly in a. The conventions are the same as in Fig. [4] 



out to be quite a bit smaller than its natural order of magnitude of 0(1/ z). In such a 
situation 0(1/ z 2 ) terms may be numerically comparable and the expansion may 

be compromised in the matching step. We therefore chose the described solution, even 
if it is lacking elegance. H The continuum extrapolation of $5(L2) is shown in Fig. [7J 

Since the computed $(£2, M, 0) may be used in the subsequent step eq. ( |3.8[ ) also 
with lattice discretizations which differ from ours we list $i(L 2 , M, 0) in Tableland Ta- 
ble |4} Starting from these numbers, the remaining computations to obtain the HQET 
parameters for a different lattice action require a very modest effort. 

6 We have of course also determined 0)5 from $s(Li). Despite the potentially large 1/m 2 contribu- 
tions, the result for £5 changes only by about 15%. 
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Figure 7: Continuum extrapolation of $5(^2) for our standard set of 9i. The central quark mass is 
shown. The conventions are the same as in Fig. [4] 





$f at (L 2 


M,0) 


z 




<j,stat 


10.4 


15.73(20) 


0.538(4) 


12.1 


17.82(22) 


0.552(4) 


13.3 


19.29(24) 


0.560(4) 



Table 3: Observables <£>f at (I/2, M, 0) for the different values of the quark mass, and for the standard 
choice of 8's. 



<ff /m) (L 2 ,M,0) 


z 








$ (l/m) 




10.4 
12.1 


0.14(7) 
0.12(6) 


0.038(14) 
0.034(13) 


0.0491(23) 
0.0430(22) 


0.0333(4) 
0.0297(4) 


0.00234(16) 
0.00208(14) 


13.3 


0.11(6) 


0.031(12) 


0.0396(21) 


0.0276(4) 


0.00193(13) 



Table 4: Observables $ 
choice of #'s. 



(l/m) 



(L2, Af, 0) for the different values of the quark mass, and for the standard 



4.3 HQET parameters 

4-3.1 Renormalization group invariant b-quark mass 



For each value of z = ML\, eq. (3.8) yields the desired HQET parameters. However, for 
future use we want to list them for M = Mb. This saves space here and somebody doing 
a computation in the future does not have to fix the right quark mass. The b-quark 



mass was already obtained in 17 from the experimental spin-averaged B s meson mass 
and rn = 0.5 fm. Here we repeat its determination using the mass of the pseudoscalar 
B s meson as experimental input. This is the natural quantity since we also start from 



15 



the effective mass defined from the /a correlator to fix the b-quark mass in the small 
volume. Moreover, we have improved some of the necessary steps through the finer 
resolution ajL<i = 1/32, and the use of tree level improvement (see App.fD]). With our 



new, improved determinations 33 using the GEVP method 22 we found P] 

r M|f at = 17.12(26) , r M^ QET = 17.38(28) . (4.3) 
These results are in perfect agreement with the ones obtained when using the large 



volume numbers E^ m , E st&t of 17 . In the following we use eq. (4.3) and the knowledge 
of L\/rQ and interpolate all results quadratically in z to = 12.30(19) in the static 
approximation and to z\, = 12.48(20) at first order in 

4-3.2 Bare parameters oji 



Our u!i determined from eq. (3.8) and simulations S4 are listed in Table [5] for our stan- 



dard 6— combinations. The errors take all sources into account (through a jackknife 
analysis incorporating all steps), but one has to be aware that there are very significant 
correlations in the parameters. We discuss these correlations in App. [C] and provide 
them in tables available on the Internet. In static approximation a small shift propor- 
tional to u\ — u\ (see the discussion at the beginning of this section) is applied. We refer 
to Appendix D of Jl7| for the case of the b-quark mass, and with similar considerations 
we have evaluated the effect on the current renormalization, ln(Z A tat ). We found 

9 = , r AM b = -0.042(20) , A ln(Z A tat ) = 0.007(5) , (4.4) 

6o = l/2, r AM h = 0.009(11) , Aln(Z A tat ) = 0.008(5) , (4.5) 

9 = 1 , r AM b = 0.150(45) , A ln(Z A tat ) = 0.010(4) . (4.6) 

Generically the bare parameters oj{ are completely non-universal and depend on 
all details of the action. However here we are working in the quenched approximation. 
In this situation the HQET action can in principle be determined independently of the 
light quarks. Thus ti>i,a/4 and 0J5 can be used for any light quark action also different 
from our specific one. 

We illustrate the cutoff-dependence of ui\ = rribare as a function of L±/a = -LiA cuto ff. 
In Fig. g we show the static bare quark mass and its l/m\, contribution in units of Li, 
using Li/a = L2/{2a) and the data of Table [5| In addition, smaller lattice spacings are 



covered by including the numbers for u)\ from eq. (3.4). The two sets of bare parameters 
u>i and Ui differ by cutoff effects: the former are determined directly for L = L\ and 
the latter after a step scaling to L = L2. Indeed these cutoff effects are visible but 



7 The reader might wonder why for roiWb we don't obtain a better precision compared to the result 



quoted in 17 . The reason is that the uncertainty which affects the b-quark mass is dominated by the 
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renormalization of the quark mass in QCD 

8 These shifts have to be added to the raw numbers computed with our simulations Si, ■ ■ ■ , S4 but 
are already included in the numbers given in eq. (4.31 and in Table [5] 
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HYP1 






HYP 2 




p 


6.4956 


6.2885 


6.0219 


6.4956 


6.2885 


6.0219 


stat 
am barc 


0.964(12) 


1.324(17) 


2.054(25) 


0.990(12) 


1.352(17) 


2.083(25) 


ln(Zf at ) 


-0.182(5) 


-0.171(5) 


-0.141(4) 


-0.118(5) 


-0.101(5) 


-0.061(4) 


(1/m) 
« m barc 


-0.315(6) 


-0.264(6) 


-0.214(8) 


-0.328(6) 


-0.273(6) 


-0.215(8) 


ln(zi 1/m) ) 


0.180(28) 


0.156(24) 


0.121(20) 


0.068(27) 


0.058(24) 


0.039(20) 


c«/« 


-0.17(7) 


-0.06(5) 


0.03(4) 


-0.61(7) 


-0.46(5) 


-0.28(4) 


^kin/a 


0.550(9) 


0.437(7) 


0.328(5) 


0.553(9) 


0.439(7) 


0.330(5) 


Wspin/ O 


0.76(5) 


0.59(4) 


0.43(3) 


0.87(6) 


0.71(5) 


0.55(4) 



Table 5: HQET parameters as a function of the bare coupling for the actions HYP1 and HYP2, for our 
choice of and for z = MhL — 12.30 (static) or z — MhL — 12.48. The three values of j3 correspond 
to L% I a — 16, 12, 8 respectively. 



not dominating. The largest part of the variation with a is due to the divergences. 
In the static case, the divergence is known to one-loop order from [18] (see table 1 of 
that paper). It is plotted for HYP1 and HYP2 using both standard (dotted lines) and 
boosted (dashed lines) perturbation theory, fixing the constant piece at the smallest 
lattice spacing. The non-perturbative results agree qualitatively with the perturbative 
approximations . 

We show the equivalent plots for ]n(Zj£ at ) and hi(Z^ m ') in Fig. § For ln(Z^ 1/m) ) 
the theoretically expected 1/a divergence is not clearly visible. It either has a rather 
small coefficient or it is masked by terms with positive powers of the lattice spacing. 
For illustration we nevertheless fit a linear behaviour to ln(Z^ m ' ) ) at the smallest three 
lattice spacings and extend this curve to larger ones. For ln(Z^ tat ) we show the leading 
logarithmic divergence, ln(Z| tat ) ~ \~n.(a / L\) g$ / {Air 2 ) + constant, in the graph. Again 
the constant is adjusted to the data point at smallest a and replacing g$ with a boosted 
coupling defines the boosted perturbation theory expression. 

Before closing this section we would like to add a remark concerning different 
matching conditions. By construction, the observables $ depend on the values of the 
angles 9{. Nevertheless, since the HQET parameters Ui are bare parameters of the 
HQET Lagrangian and fields, they are 9 independent up to truncation corrections of 
order (l/m^) 11 . As one can see from eq. (3.1 ), this implies that the ^-dependence of $ is 
absorbed by the ones of r\ and <j). In practice this means that a parameter u computed 
in the static approximation, such as or Zf^ exhibits a small 0-dependence, but 

once the l/m^ corrections are added, this dependence has to be further suppressed to 
order (l/m^,) 2 . We checked that we obtain this behaviour in our numerical simulations. 
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Figure 8: Static contribution (left) and 1/mb correction (right) to the heavy bare quark mass. In the 
static case, we also show a comparison with standard and boosted perturbation theory [l8], as described 
in the text. Data for u)\ from Table [5] (large volume simulations S4) are represented by circles (the three 
points on the left), while Q\ obtained from the small volume simulations S2 are represented by squares. 
We use the colour blue for HYP2 and red for HYP1 (also slightly shifted to the right). The results are 
shown for the central mass and for the standard set of 8i. 



-0.05 
-0.1 

+^ 

^< 

N -0.15 
-0.2 
-0.25 



0.6 

X< 0.2 . , t 

^ ...4— -t— t--f- 

M * 

-0.2 



10 

Li/a 



15 



10 

Li/a 



15 



Figure 9: Same as figure Fig. H for \n(Z A ). 
It is illustrated for the case of 0J4 in Table [6j 





L 2 /a 


(61,62) = (0,0.5) 


(6 1 ,6 2 ) = (0.5,1) 


(01,0 2 ) = (0,1) 






0.0098(3) 


0.0297(4) 


0.0395(7) 




16 


0.566(18) 


0.565(7) 


0.566(10) 


Wkin/a 


12 


0.450(15) 


0.449(6) 


0.450(8) 


Wkin/ a 


8 


0.339(11) 


0.338(4) 


0.338(6) 



dependence of $4 m '(L2, M, 0) and the corresponding HQET parameter u>i(M, a) 



Table 6: 

widn(M, a), for 2 = 12.1 and the action HYP 2 
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5 Outlook 



Our non-perturbatively computed HQET parameters show a qualitative agreement with 
perturbation theory as far as the leading divergence in each of the parameters is con- 
cerned (see figures 8|9 ), but this does not extend to the quantitative level needed for 



precision flavour physics. 

Fortunately, in the quenched approximation, we now have the full set of parameters 
for HQET spectrum calculations as well as for (zero space-momentum) matrix element 
of Aq including the terms of order 1/mb - all of them are known non-perturbatively. 
A detailed test of the 1/m^ expansion is thus possible. In companion papers we are 
carrying this out for the examples of the B s decay constant and for some mass splittings 
in the B s system; some preliminary results are described in 1 33 1 . The parameters do 
of course depend on our choice of discretization, but starting from Tables [3] and |4j the 
effort to compute them for a different regularization is very modest. Note that, as long 
as one remains in the quenched approximation, the parameters for the HQET action do 
not need to be recomputed if one uses just a light quark action differing from ours. 

The small volume simulations for a determination of the parameters for Nf = 2 
are already far advanced, see [7] for a recent account. Therefore, the situation will soon 
be similar with 2 flavours of dynamical fermions. 
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lar H. Simma for very useful comments on the manuscript. This work is supported by the 
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financial support from the MICINN grant FPA2006-05807, the Comunidad Autnoma 
de Madrid programme HEPHACOS P-ESP-00346, and participates in the Consolider- 
Ingenio 2010 CPAN (CSD2007-00042). 

A Choice of observables 



Our observables are built from Schrodinger functional |34 35 correlation functions de- 
fined exactly as in [IT] . We use the renormalized boundary-to-bulk correlation function 
/a(xo,#) m * ne pseudoscalar channel as well as the boundary-to-boundary correlators 
fx, k\, in pseudoscalar and vector channel respectively. For their interpretation, we note 
in particular the role of 9. Quark fields ip are periodic in space up to a phase e td , while 
ijj are periodic up to e~ lS . This shifts the lowest momentum of a free quark from p = 
to p = (1, 1, 1)6 /L. Intermediate states to the above correlation functions have vanish- 
ing momentum, due to a projection a 3 ^ x and the fact that fields ipiFfa are periodic 
(with a vanishing phase). However, in the free theory, which is a relevant starting point 
for the interpretation in small volume, the quark carries momentum p = (1, 1, 1)8 /L, 
compensated by the opposite momentum of the anti-quark. 

We form an effective energy T p = — 9o(ln(— /a(^o> $o)) which ap- 

(x =T/2,T=L) 

proaches the heavy quark mass as it becomes very large. The logarithm of the ratio 
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i?A = 1h(/a(T/2, ^i)//a(T/2, 02)\t=l 1S easily seen to have a sensitivity to the coeffi- 
cient W3 which is approximately proportional to 9 2 — 9\\ one just has to note that the 



covariant derivatives in eq. (2.5) acting on the quark fields are proportional to their 
momentum. This free theory argument is valid qualitatively in small volume. In the 
same way one sees that the combination R x = \ (ln(/i(0i)/;i(0i) 3 ) - Mh{e 2 )k 1 {e 2 f)) 
has a sensitivity proportional to Q\ — Q\ to W4, while W5 does not contribute due to spin 
symmetry. 

These qualitative considerations combined with some numerical experiments lead 
us to introduce 

* = (Lr p , In ,R A , Rl , 3 -ln (£) )* . (A.l) 



B Explicit form of step scaling functions 

From our definitions we find immediately 

$(L 2 ,M,0) = (B.l) 
lim (7/(^,0) + 0(L 2 ,a)^ 1 (Li,a)[$(L 1 ,M,O) -77(^,0)]} , 

a/Li->0 



where 



, (T nU - 1(T , (D [B(L 2 , a) - DB(Li,a)] a)\ 

D = diag(L 2 /Li, 1) (B.3) 



in terms of the matrices introduced in eq. (3.5). The form 

$;(L 2 ,M,0) = A*i(£i,M,0)+ lim %{L u a) (B.4) 

a/Li— >0 



+ lim 0)- 7]j{L x , 0)], i = 1,2, 



J=3 



5 



$i(L 2 ,M,0) = 7/^2,0)+ lim V^^.M.ajftfLi.M.Oj-^lLi.O)], 

z = 3...5. (B.5) 

then introduces step scaling functions 

Sj(Li,a) = r)i(L 2 , a) - Di^Lx, a) , i = 1,2 , (B.6) 
E(L, a ) = (° W^a)-fl^,.)M-'(Ll..)y (B.7) 
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By choosing only L\ and a as arguments we have assumed that L2/L1 = s is fixed 
(typically to s = 2) which means D = diag(s, 1). 

The continuum limit of each element of the step scaling functions exists and the 
above split into two pieces is suggested by the fact that the limit a — > of r]i(L, M, a) 
may be performed for i > 3, while for i < 2 there is an additive renormalization which 
only cancels in eq. (B.6). Splitting r] accordingly, r/ = rf + r/ b with rif >3 = , r/^ <2 = 0, 
we can also unify the step from L\ to L2 to (note S(L,a)j>3 = 0) 

$(L 2 ,M,0) = V \L 2 ,0)+ lim £(Li,a)[$(Li,M,0) -r/ b (Li,0)] +%{L l ,a). (B.8) 

a/Li^O 

Let us now list explicit expressions for the various matrices. The expansion of the 
observables in HQET follows directly from the expansions detailed in section 2.3 of |17| . 
The inhomogeneous part is given by the static quantities 



where 



r stat (L) 

Ca(L) 
Rf at (L) 
i?f at (L) 



(Lr stat , Ca, #A a \ Rf at , o) 1 , 

do + d* 



(x =T/2,T=L) 



111 

In 
In 



Vf! tat (0o) . 

/ A tat (^o^i) 

^stat 
Jl 

J 



(x =T/2,T=L) 
(x =T/2,T=L) 



stat 1 



(B.9) 

(B.10) 
(B.ll) 

(B.12) 
(B.13) 



T=L/2 



Note that / A tat does not contain the improvement term proportional to c^ at . The 
improvement term is included as explained in section |3j From the definition eq. (B.6), 
one has 



Si 



L 2 (r stat (L 2 )-r stat (Li)), 
Ca(l 2 ) - Ca(Li)) • 



(B.14) 
(B.15) 



The matrices A, B which make up (f> as in eq. (3.5) are 



A 




B 



LTsa LF kin LF s P' m \ 
Psa ^ kin * spin / 



(B.16) 
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where 



I*(L) 
Tsa{L) 

P Y (£) 
Psa{L) 
R\{L) 

Rsa(L) 
R\(L) 



fi(x ,e ) 1 fl(9 ) 



.f A tat (*oA) 2f! tat (0 ] 
fl(Oo) 



(x =T/2,T=L) 
(x =T/2,T=L) 
{x =T/2,T=L) 



A stat (#o) 

/Ja^oA) 



T=L/2 



/ A tat (x O ,0o) ( XQ =T/2,T=L) 

/IMi) /I(x ,g 2 ) 

/ A tat (^0^l) / A tat (^0^2) 
/l tat (^0^l) / A tat (^0,^2) 



(o;o=T/2,T=L) 
(a;o=T/2,r=L) 



/f(fll) 



A 7 



T=L/2 



In order to avoid unnecessary repetitions, we have used the shortcut 

y G {kin, spin} . 
in the previous formulae. From our definitions we obtain 



A-' L = 



Kl2 = 



I/Rsa K 12 K 13 
1/Rf n 
l/pf n 

-Rt n 



pkin 

-K A 



R\ m R8A 



, K 13 = 



and 



[5(L 2 ,o)-DB(Li I op- 1 (Li, a ) 



Pi ^<5A 



S13 S14 S15 

S23 S24 S25 

S33 S34 £35 

s 44 

£55 
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with 



S13 


= L 2 (Tsa(L2. 


) - r SA (L x ))/R SA (L x ) , 


S14 


= L2{T 5A {L 2 


) - r M (Li)) K 12 {L X ) + L 2 (r kin (L 2 ) - r kin (L x ))/ R x in (L%) , 


S15 


= L 2 (Tsa(L 2 


) - r SA (L x )) K 13 (ii) + L 2 (r spin (L 2 ) - r spin (L 1 ))/ P f n (L x 


^23 


= (psa(L 2 ) - 




^24 


= (psa(L2) - 


m(^i)) ^12(^1) + (* kin (£ 2 ) - ^(loj/^cli) , 


^25 


= (psa(L2) - 


Psa(Li)) K X3 (L x ) + (* spin (L 2 ) - ^ spin (L 1 ))/ / of n (L 1 ) , 



(B.30) 

and 



£33 


= R SA (L 2 )/R SA {L X ), 


S34 


= R$a{L 2 )K\ 2 (L\) + i?i in (L 2 )/ J R kin (L 1 ) , 


£35 


= R SA (L 2 )K X3 (L X ) + if/ n (L 2 )/pf n (Li) , 


£44 




£55 


= p? m (L 2 )/Rf n (Li). 



(B.31) 



In static approximation the only non-vanishing elements are 



13 — ^13 > ^23 — ^23 ■ 



C Covariance matrices 



When our parameters are used in a subsequent computation, one has to note that their 
errors are correlated. We therefore supply covariance sub-matrices as far as they appear 
necessary in practice. They are still too large to be printed (and typed in). We therefore 
supply them for download together with the paper under http : / / arxiv . org/ and at 
http://www-zeuthen.desy.de/alpha/public_tables/. All parameters are given in 
lattice units a = 1, and for our standard choice of 6 angles. As explained in the text, we 
give the HQET parameters and their covariance matrices for c^ at fixed to its one-loop 
value pO] and fixed to 0. 

We normalize the covariance matrix between two quantities O a and O/3 in the 
following way 

r - Ga P 



(Tafl = {N iaA -l){{O a Op)-{O a ){Op)), (C.l) 

where ( ) represents the average over = 100 jackknife samples. The covariance 

matrix is such that Cp a = C a /3 and C aa = 1 (no summation). 
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A first subset of parameters is (wf tat , u[ , 014), necessary for the computation of the 
quark mass. In order to take into account the different kinds of correlations (correlation 
between different values of the lattice spacing, or between different values of the heavy 
quark mass), the indices a and (3 take the values 

a, /3 = I + 3(i - 1) + 9(j - 1) + 27(k - 1) . (C.2) 

Here I = 1,2,3 represents the parameters wf tat , oj^ and W4 = u^ m ^ respectively, 
i = 1,2,3 labels the three go values, j = 1, 2, 3 numbers the three values of the quark 
mass and k = 1(2) for the HYP1(2) static actions. This leads to a 54 x 54 matrix C^g, 
but for symmetry reasons, we give only the 1431 components with a < f3. 

The second subset is the one relevant for a computation of the pseudoscalar decay 
constant. There the indices a and j3 take the value 

a,0 = J + 5(*-l) + 15(fc-l). (C.3) 

Here I = 1 represents w| tat and I = 2 ... 5 represent co^ , . . . ,co^ m \ while i, k have 
the same meaning as before. Thus we obtain a 30 x 30 matrix C*g, and give its 435 
relevant components. The heavy quark mass is fixed by matching to the spin averaged 
B s -meson mass. This matching is done at the static order for u; stat and including the 
1/mb terms for aj( 1//m ). 



D Tree level improvement 

Perturbative improvement of the observables [36] has been proven to be effective in 
reducing cutoff effects for example in the case of the step scaling function for the cou- 
pling [37j[38], where it has been pushed to the two-loop order. The idea is to remove 
all the 0((j^) n ln(-j) m ) at a given order in the g 2 expansion of the renormalized lattice 
observable. At tree level, as done here, this produces classically perfect observables in 
the sense of 39 , by removing all the 0((^) n ) effects. 



The tree-level improved observable is either defined as 



O impr {a/L) 



Q(a/L) 
1 + 5 (a/L) 



8 (a/L) 



O tree (a/L) -O tree (0) 
O tree (0) 



(D.l) 



or as 



O 



impr 



(a/L) = O (a/L) - 5 (a/L) , 5 (a/L) = O tiee (a/L) - O tree (0) . 



(D.2) 



Here O is the non-perturbative observable evaluated by Monte Carlo, Oimpr is the tree- 
level improved one and O tree is the same observable evaluated for go = 0. Generically 



we use eq. (D.l), but when O tree (0) vanishes eq. (D.2) is appropriate. We also choose 



the latter for the step scaling functions S since it is more natural given their form 



eq. (B.14), eq. (B.15). We apply this improvement to all QCD and HQET observables 



used for the matching in small volume and to all HQET step scaling functions. 
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At tree level there are some useful relations among the Schrodinger Functional 
correlation functions, which can be used in the perturbative computation |40| : 



ky(x Q ) 

h 



3 
fei 



/p(^o) - o/a^o) , 



:fp(xo) - o/aC^o) 



(D.3) 
(D.4) 



XQ=T—a 



Since in our QCD simulations L is fixed through g 2 (L), each value of z determines 
a value of M. In the corresponding tree level computation the improved subtracted 
quark mass^jmq = m q (l + 6 m am q ) enters. At the considered order we could set m q to 
m(n) at an arbitrary scale /i in an arbitrary scheme. We decided to choose m q = m*, 
where m* = m un i v (m Vr ) and fn un i v (fi) is the running mass computed from the RGI mass 
(but) using only the universal part of the f3 (i.e. the coefficients &o an d b\) and r (i.e. 
the coefficient do) functions. This choice is scheme independent and it is based on the 
expectation that cutoff effects are dominated by scales around \x = m*. The explicit 
relation between z and m*, 



A SF L 



2 &i/(26g) ^1-61/(60^0) 



-2b /d 



with 



(D.5) 



is then solved for m*. In our non-perturbative computation, the $^ CD have been 
computed for z = ML = 10.4, 12.1 and 13.3 with g 2 (L) = u\ (see the begin- 
ning of Sect. [4]) and correspondingly Asp£ = 0.195(16). Our tables, provided under 
http://www-zeuthen.desy.de/alpha/public_tables/, are for these specific values. 

As an example, the non-perturbative $^ CD with (right) and without (left) tree 
level improvement of the observable are compared for (9\, O2) = (0.5, 1) in Fig. 10 



0.3 r 



0.29 



8 0.28 
0.27 



0.26 




0.3r 
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O 0.28 
0.27 



0.26 
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x 10"'' 
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Figure 10: Effect of tree level improvement on $J CD . The three masses are shown, with the same 
conventions as in fig [3] 



Here 6 m is an improvement coefficient 41 
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Turning now to HQET, the correlation functions are denned in [17] and useful 
relations at tree level are 

Y ) " l) /I tat (^o) , (D.6) 

( T ~ a ) f„„J a0 \ 1 \ ,stat 





= 6 ^ 

a z 


^cos 


11 


= kf n 


= 6 


Jx 


= . 





cos ( -J - 11 /f at , (D.7) 

(D.8) 

The last relation holds when the magnetic background field vanishes, which is the case 
in our application. 

It is instructive to consider now the step scaling function £44, which has also been 
discussed in |17^, but without tree level improvement. Explicitly 

E u (u, a/L) = [RY*(2L)/Rf»(L)]- g2{L)=u , (D.9) 

with 

Rf n (L,ex,e 2 ) = /j^iV/f^i) - ft m {e 2 )/fr\e 2 ) , t = l/2. (d.io) 

Given the tree level relations among correlation functions, the expression for the step 
scaling function reads 

L_ 1 (cos (&) -cos 

E&>/L) = -f— ^ • ) ^ )^ , (D.ll) 

4-1 (cos cos (^ 

with S44 CC (0) = 0.5 . The first factor is responsible for cutoff effects linear in a, which 
appear because we have not taken care to 0(a) improve the theory at order l/m^. The 
correction 5 is then large and tree-level improvement leads to a strong reduction of the 



cutoff effects in the non-perturbative S illustrated in Fig. [11 



For future applications, this discussion suggests to eliminate linear (tree level) a/L 
effects such as those in Fig. 11 from the start. This is achieved by implementing 0(a)- 
improvement at tree- level in the order 1/mb terms of the action and the effective fields. 
It just means to add terms 

5S = -a 3 ^^O kin (0,x) + ^ kin (T,x)], (D.12) 

X _ 

4 3) = -i^boTsD 2 ^) , (D.13) 

to action and currents. The normalization is chosen such that a coefficient will 
implement tree- level improvement, 

^ ET (-) = ^[^rw+i^f'+i^P+^a, (D.i4) 

10 There it is called £ kin . 
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and similarly for the action. The effect of these terms is that the kinetic operator D 2 
is inserted only with weight 1/2 on the initial and final time-slice of the correlation 
functions such as eq. ( |2.12[ ) or f%*, ff n , jfej* n . After the Wick-contraction (in a given 



gauge background field) this corresponds to a standard discretized representation of the 
integral over the time position of the insertion of D 2 in the static propagator. The spin 
operator O sp i n can be treated in complete analogy, but due to eq. (D.8) this has no 
effect at tree-level. 
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r ^0.55 
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Figure 11: Effect of tree level improvement on E44 for (61,82) 
same as in fie Hl 
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(0.5, 1). The conventions are the 



We provide the tree-level improvement coefficients under 
http : //www-zeuthen . desy . de/ alpha/public_tables/. 
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